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Abstract 

Radiation of X-ray bursts and of accretion shocks in weakly magnetized 
neutron stars in low-mass X-ray binaries is produced in plane-parallel at- 
mospheres dominated by electron scattering. We first discuss polarization 
produced by single (non-magnetic) Compton scattering, in particular the 
depolarizing effect of high electron temperature, and then the polarization 
due to multiply electron scattering in a slab. We further predict the X-ray 
pulse profiles and polarization properties of nuclear- and accretion-powered 
millisecond pulsars. We introduce a relativistic rotation vector model, which 
includes the effect of rotation of polarization plane due to the rapid motion 
of the hot spot as well as the light bending. Future observations of the X-ray 
polarization will provide a valuable tool to test the geometry of the emission 
region in pulsars and its physical characteristics. 



1.1 Introduction 

Polarization has proven to be a valuable tool in determining the geometry 
of the emission region in radio pulsars |l|). For the X-ray pulsars, the data 
are not yet available, but their interpretation in any case is not going to be 
easy, because of the strong magnetic field effects on the radiation transport. 
Discovery of millisecond coherent pulsations during X-ray bursts in nearly 
20 low-mass X-ray binaries (so called nuclear-powered millisecond pulsars, 
NMSP, see (0)) and in the persistent emission of at least eight sources 



(accretion-powered millisecond pulsars, AMSP, see (|25l ; 1 14 : Il5l)) opens a 
completely new range of possibilities. The emission in these cases is pro- 
duced at the surface of a rapidly spinning, weakly magnetized neutron star. 
Thus the magnetic field does not affect the radiation transport and much 
more reliable predictions for the radiation pattern from the surface can be 
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obtained. The observed pulse profiles and polarization are affected not only 
by general, but also special relativistic effects (0; 0; 0) • The pulse profile 
alone, however, does not allow to determine uniquely the pulsar geometry, 
while the phase dependence of the polarization angle (PA) is a powerful 
tool to distinguish between the models and to put strong constraints on the 
geometry. In this review, we first discuss the physics of polarization pro- 
duced by Compton scattering. Then we apply these results to predict the 
polarization properties of NMSP and AMSP. 



1.2 Polarization properties of Compton scattering 

1.2.1 Polarization in single Compton scattering 

Radiation Thomson scattered once by cold electrons becomes linearly po- 
larized with the PD depending strongly on the scattering angle ©: P(O) = 
sin 2 0/(1 + cos 2 0). With increasing electron temperature the PD drops 
because of random aberrations and corresponding random rotations of the 
polarization plane due to the thermal motions of the electrons (12; 1^) (see 
Fig. II. ip . For isotropic relativistic electrons of arbitrary energies, the exact 
analytical expressions for the scattering redistribution matrix for Stokes pa- 
rameters was derived in (fill ). For T c = 50 keV, polarization is smaller by 
50% and for 100 keV it drops by a factor of two compared to the Thomson 
scattering case. As electrons in accreting black holes and neutron stars often 
reach such temperatures, the depolarization effect has to be accounted for. 
For isotropic power-law distribution of relativistic electrons, the net scat- 
tered polarization is zero if the incident polarization is unpolarized 




Fig. 1.1. PD at 1 keV of 0.1 keV blackbody photons scattered by the hot isotropic 
electrons of various temperatures T c as a function of the scattering angle (adapted 
from (fl3h ). For the dependence of the PD on photon energy and for cases of initially 
polarized radiation, see (QT 
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1.2.2 Scattering in optically thick atmosphere 

In a classical problem of radiative transfer in a plane-parallel semi-infinite 
atmosphere with opacity dominated by Thomson scattering (0; I21I ) , the PD 
reaches the maximum of ~ 11.7% at \i = (where \i is the cosine of the 
angle from the normal) and decreases to zero at u = 1 due to the symmetry: 

Pra 1 ~ M 11.71%. (1.1) 

1 + 3.582/j V ; 

The dominant direction of the electric vector oscillations is parallel to the 
slab plane. 



1.2.3 Comptonization in optically thin atmosphere 

The detailed study of the transport of polarized radiation in plane-parallel 
atmospheres with absorption and Thomson scattering included was consid- 
ered in (d; 0) . The absorption causes the rotation of the polarization plane 
by 90° (Nagirner effect, did)) even in optically thick atmospheres. The same 
happens when the scattering optical depth becomes small tt < 1 (|22j; |6|; 1 24 ) . 
We are often interested in polarization of photons undergoing a certain num- 
ber of scattering (because photons gain energy every time they are scattered 
by hot electrons). The dependence of intensity and PD on angle relative to 
the normal for various scattering orders is shown in Fig. 11.21 for tt = 1. 
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Fig. 1.2. Angular dependence of the intensity and polarization of the radiation 
escaping from a slab of Tr = 1 for different scattering orders n. Positive P corre- 
sponds to the plane of polarization parallel to the slab normal. The unpolarized 
seed photons with =const were injected from the slab bottom. The dashed 
curves correspond to the classical results of Chandrasekhar-Sobolev for tt = oo. 
From dU 
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Usually polarization was considered in Thomson scattering approxima- 
tion, while, as we have seen in £ jl.2,ll the polarization does depend on the 
electron temperature. Calculations of the PD in a slab using fully relativis- 
tic kernel from ([ill ) is shown in Fig. 11.31 One sees that radiation scattered 
once is polarized perpendicular to the slab normal (as for optically thick at- 
mosphere), while at higher scattering orders polarization changes the sign. 
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Fig. 1.3. The flux and the PD emergent from a slab of tt = 0.5 and electron tem- 
perature kT e = 55 keV at inclination of 60°. The incident photons are blackbody of 
temperature 10 eV. Different scattering orders are marked at corresponding curves. 
The thick solid curves correspond to the total radiation. The zeroth order radiation 
consists of unscattered blackbody and a polarized Compton reflection component 
at high energies (computed using Green's matrix from |l9|)). Adapted from (flU ). 



1.3 Polarization of accreting neutron stars 

1.3.1 Relativistic rotating vector model 

When considering polarization of radiation from coherently oscillating sources, 
we can assume that the emission originates in one or two (in the case of 
AMSP) spots at the neutron star surface. In AMSP the hard X-rays are pro- 
duced by thermal Comptonization in a plane-parallel slab (accretion shock) 
of optical depth of order unity dl7|; Il5l ) . The emission during the X-ray bursts 
is produced in a semi-infinite electron scattering dominated atmosphere. It 
is natural to assume that the radiation pattern has azimuthal symmetry in 
the comoving frame of a spot. Once we have the Stokes parameters in spot 
frame, we can transform them to the observer frame. First, we make the 
Lorentz transformation to the non-rotating frame accounting for Doppler 
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boosting and relativistic aberration and then follow photon trajectories to 
the observer at infinity in Schwarzschild space-time. Deviations from the 
Schwarzschild metric and from sphericity of the star due to the stellar ro- 
tation have a small effect and are usually neglected. For pulsar rotational 
frequencies of v > 400 Hz, we also need to account for time delays which in 
the extreme cases can reach about 5-10 per cent of the pulsar period. 

Lorentz transformation and gravitational light bending do not change the 
PD and the observed value corresponds to the polar angle a' at which a 
photon is emitted in the spot comoving frame. Lorentz transformation gives 
cos a' = 5 cos a, where 5 is the Doppler factor and a is the angle in the 
non-rotating frame, which is related to the position angle ip by the usual 
gravitational bending formula (see Fig. 11.41 and fl7l ; I24I: llq) for details). 
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Fig. 1.4. Geometry of the emitting regions for AMSP. 



For a slowly rotating star, the polarization vector lies in the plane formed 
by the spot normal and the line of sight. The PA (measured from the 
projection of the spin axis on the plane of the sky in the counter-clockwise 
direction) is given by 

sin# s'm<j) 

tanxo = — — 2 — ~a 1' C 1 - 2 ) 

sm 1 cos — cos 1 sm cos <p 



as in the rotating vector model of Radhakrishnan & Cooke (20j) . Here cf> is the 
pulsar phase, 6 is the spot colatitude, and i is the observer's inclination (see 
Fig. II. 4p . The effect of the stellar spin on rotation of the polarization vector 
is discussed in (|4j; |5j) . Viironen & Poutanen (|24j) presented the correction to 
the expression for the PA which includes light bending and aberration: 

. cos i sin 6 — sin i cos 9 cos <f> 

tanxc = Peqcosasinfl- : — — — (1.3) 

sm a sm ip + p eq sm o sm 1 sm <p 
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Fig. 1.5. Light curves, PD and PA expected from an X-ray burst happening at 
various colatitudes 8. Semi-infinite electron scattering atmosphere is considered. 
The neutron star mass is M = 1.4M Q , radius R = 10.3 km, rotational frequency 
v = 400 Hz, and the observer inclination i — 60°. From (1241). 



where /3 eq is the equatorial stellar velocity in units of speed of light. The 
total polarization angle for each spot is then x = Xo + Xc- 



1.3.2 Polarization properties of NMSP 

The X-ray bursts often show coherent oscillations (23). The energy dissipa- 
tion takes place deep in the atmosphere, and thus we can assume that optical 
depth is infinite. At sub-Eddington luminosities the atmosphere is pinned 
down to the neutron star surface, and thus a plane-parallel approximation 
is valid. At effective temperatures of about 2 keV, most of the material is 
ionized and scattering dominates the opacity. The radiation escaping from 
the surface then can be described by Chandrasekhar-Sobolev formulae from 
£ 11.2.21 The predicted pulse profiles and behavior of PD and PA are shown in 
Fig. ll.5l for a small bright spot. The polarization is increasing with the spot 
colatitude 6, reaching the maximum of ~ 12% close to the eclipses. The PA 
varies around 90° as the electric vector is predominantly perpendicular to 
the meridional plane and its variability amplitude grows with 9. The larger 
is the spot, the smaller is the polarization. 
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Fig. 1.6. Normalized pulse profiles, PD and PA from two antipodal point-like spots 
for various inclinations and magnetic inclinations. The intrinsic radiation is scat- 
tered 7 times in a slab of Thomson optical depth tt = 1 . Dashed curves correspond 
to the primary spot, dotted curves to the antipodal spot, and solid curves to the 
total signal. The stellar parameters are M — 1.4M Q , R = 10.3 km, v = 300 Hz. 
From j2J). 



1.3.3 Polarization properties of accreting millisecond pulsars 

The spectra of accreting millisecond pulsars can be represented as a sum 



of a black body-like emission and a Comptonized tail (see review (|15l ) and 
references therein) . The black body photons play also a role of seed photons 
for Comptonization. This interpretations is supported by the shape of the 
pulse profile above 7 keV which is consistent with being produced by the 
hotspot with the angular pattern characteristic for an optically thin slab 
0). Because below a few keV contribution of the black body is large 
the expected polarization degree is small. Scattering in the hot electron 
slab modifies significantly the angular distribution of radiation and produces 
polarization signal as described in § 11.2.31 The polarization degree is a strong 
function of the scattering order for small n (see Fig. 1 1 . 3j) and therefore of 
the photon energy also. After a few scatterings polarization saturates. On 
Fig. 11.61 we show the prediction for the polarization signal for n = 7 (in 
Thomson approximation), which would correspond to about 10 keV for a 
typical electron temperature of 50 keV and the seed photon temperature of 
0.6 keV. We choose the slab optical depth tt = 1 which is consistent with 
the spectra of SAX J1808.4-3658 |l3;Q). We see that exchanging i and 6 
does not affect the pulse shape and PD, but the PA changes dramatically. 
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1.4 Conclusions 

Polarization properties of non-magnetic Compton scattering, which is the 
main emission mechanism in weakly magnetized accreting neutron stars, 
have been studied in details. Light bending in the gravitational field of the 
neutron star and the rotation of the polarization plane because of relativistic 
effects are also well understood. Thus the theory of the polarization from 
millisecond X-ray pulsars is rather well developed by now and is ready to 
be applied to the data. The launch of an X-ray polarimeter will hopefully 
start a new era in studies of these sources and will provide a powerful tool 
for understanding the geometry and the physics of pulsars. 
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